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We apply the fermion functional renormalization group method to determine the pairing sym- 
metry and pairing mechanism of the FeAs-Based materials. Within a five band model with pure 
repulsive interactions, we find an electronic-driven superconducting pairing instability. For the dop- 
ing and interaction parameters we have examined, extended s-wave, whose order parameter takes 
on opposite sign on the electron and hole pockets, is always the most favorable pairing symmetry. 
The pairing mechanism is the inter Fermi surface Josephson scattering generated by the antiferro- 
magnetic correlation. 



The copper oxide high temperature superconductors 
(the cuprates) have stimulated a significant part of con- 
densed matter physics development in the past twenty 
years. Two years ago, a new non-copper superconduct- 
ing compound, the iron pnictides, was discovered [1]. In 
the last few months, by various element substitution, 
the superconducting transition temperature of this new 
class of superconductor has been raised to 55K[2, 3]. 
This has stimulated a flurry of interests in their ma- 
terial and physical properties. Examples of iron pnic- 
tides include LnFeAsOi_a;Fa,[l-3], Lai_a,Sra,FeAsO[4], 
and Bai_a;Ka;Fe2As2[5]. 

Unlike the cuprates, the stoichiometric parent com- 
pounds of the iron pnicides are metallic (or semi- 
metallic) [6] antiferromagnets[7] rather than antiferro- 
magnetic Mott insulators. For example, recent angle- 
resolved photoemission spectroscopy on single crystal 
BaFe2 As2 has revealed residual Fermi surface in the anti- 
ferromagnetic state [8]. On the other hand, similar to the 
cuprates, the antiferromagnetic order is quickly replaced 
by superconductivity as a function of doping [9], suggest- 
ing an intimate relation between antiferromagnetism and 
superconductivity. Currently it is generally felt that pro- 
gresses in understanding this new class of high Tc super- 
conductor will also have impacts on the two decade old 
problem - the cuprates. 

If one reviews the physics of the cuprates, it is easy 
to spot the prominent role played by the pairing sym- 
metry. At the present time early experiments on the 
iron pnictides have produced conflicting claims on this 
issue [10, 11]. As far as theoretical proposals are con- 
cerned, pairing symmetries spread all over the map 
[12, 13]. 

The purpose of this paper is to determine the pairing 
symmetry and pairing mechanism theoretically. How- 
ever, because of the intermediate coupling nature of these 
material [14], it is not obvious where to start. For exam- 
ple, if one takes a model of the bandstructure and adds 
repulsive local interaction to it, it is easy to show that at 
the mean-field level there is no superconducting instabil- 



ity in reasonable parameter range. In addition, unlike the 
cuprates, there is no strong coupling (the t^ /U ) expan- 
sion that can lead to antiferromagnetic exchange which 
in turn can produce pairing (at least at the mean-field 
level). In the process of searching for a suitable calcu- 
lational scheme we come across the functional renormal- 
ization group (FRG) method[15]. The numerical version 
of this method was implemented in Ref. [16, 17] in the 
study of the cuprates. Interestingly it successfully repro- 
duced both the antiferromagnetic and the d-wave pairing 
tendencies. The advantage of this method lies in its abil- 
ity to generate effective interactions which are otherwise 
absent at the bare level. In addition, by monitoring the 
growth of the renormalized interaction, it is possible to 
pin down the cause of different types of order. Finally 
in view of the fact that the electron-electron correlation 
of the iron pnictides is weaker than that of cuprates [14], 
we felt that the FRG method will have a good chance to 
succeed. 

To model the band structure, we take the most com- 
plete tight binding model we can find[18], then we add 
the Hubbard-like and Hunds-like local interaction to it. 
Like most works in the literature we focus on the As-Fe- 
As tri-layers, and view the arsenics as mediating hopping 
between the five iron 3d orbitals. The resulting tight- 
binding Hamiltonian reads 



^o = EE^lk.^'^''(kh 



6ks 



(1) 



k,s a, 6=1 



where c^ks annihilates a spin s electron in orbital a 
and momentum k.The parameters used in constructing 
Kahi}^) can be found in Ref. [18]. However, in our cal- 
culation a gauge has been chosen so that all elements of 
Kahi)^) are real. Doping is controlled by adding a chem- 
ical potential term to Eq. (1). For sufficiently large elec- 
tron doping, there are two different hole pockets centered 
at k = (0, 0) and two electron pockets centered at (tt, 0) 
and (0,7r). For hole doping, or small electron doping, an 
extra hole pocket is present at k = (tt, tt). These features 
are consistent with other DFT calculations [13, 19]. 



We describe the electron-electron interaction by 



^int =^{UlYl ^i,«,T^i,«,i + ^2 XI rii^arii^p ^Jh[^Y1 4as4f3s'^ias'^iPs + {4a}4ai^i(3iCi(3} + h.C.)\ }. (2) 



a<f3 



a<f3 s,s^ 

r 



Here i labels the sites of a square lattice, s, s^ =t, |, and 
^2, a = ^i,Q;,T + ^i,a,i ^^ ^^^ number operator associated 
with orbital a. This Hamiltonian includes the intra and 
inter orbital Hubbard Ui and /72, the Hund's interaction 
Jh, and the inter orbital pair hopping. The total Hamil- 
tonian H = Ho -\- Hint is the starting point of our study. 
The bare interaction parameters we use throughout the 
rest of the paper are Ui = 4.0, U2 = 2.0, and Jh = 0.7 
eV. The energy scale 4eV is taken from Ref. [14], and 
we have chosen the parameters so that it approximately 
satisfies the relation Ui = U2 -\- 2Jh [20]. 

We have checked that mean field theory done on the 
bare Hamiltonian, Ho-\-Hirit, has no superconducting in- 
stability for realistic interaction parameters Ui^2 > Jh- 
In following we shall show that as the high energy elec- 
tronic excitations are recursively integrated out an ef- 
fective interaction that drives extended s-wave pairing is 
generated. 

Many of the technical details of our FRG can be found 
in Ref. [17]. At the end of the present paper we dis- 
cuss few technical challenges in generalizing the method 
to the present multi-band and multi-Fermi surface situ- 
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FIG. 1: (a) A schematic representation of the gap functions 
on different Fermi surfaces. The Fermi surfaces are shown 
in black dashed lines. The width of the region between the 
solid line and dashed line indicates the magnitude of the gap 
function, and the sign is represented by the colors, (b) Two 
(red and blue arrows) typical inter Fermi surface pair transfer 
(Josephson) processes that drive pairing, (c) The five Feyn- 
man diagrams contributing to the renormalization of the ver- 
tex function. 



ation. In brief, we divide the Brillouin zone (BZ) into 
N patches, and at each renormalization iteration we sum 
over the five one- loop Feynman diagrams in Fig. (l)(c) 
in computing the renormalized 4-point vertex function 
F(ki,a;k2, 6; ka, c; k4, d). Here a, 6, c, (i = 0, . . . ,4 labels 
the five different bands. The momenta ki,..,k3 are on 
the Fermi surfaces. (Due to momentum conservation the 
fourth momentum k4 is determined.) The convention is 
chosen so that the spin associated with ka (k4) is the 
same as that associated with ki (k2). Like in Ref. [17], 
approximations (such as ignoring the frequency depen- 
dence of the vertex function, and projecting the external 
momenta onto the fermi surfaces) are made in our calcu- 
lation due to practical limitations. 

From the renormalized vertex function we extract the 
effective interaction in the Cooper channel as follows 



K1(k,a;p,5) 

= V^(k, a; — k, a; p, 6; 



-p, b) ± V(-k, a; k, a; p, b; -p, b). 



Here the upper/lower sign is for singlet/triplet pairing, 
respectively. After the BZ discretization k and p only 
take on a finite number of values. Therefore we can treat 
Vff as matrices with (k, a) and (p, b) as indices. The 
eigenvalues of this matrix are the effective interaction 
strength in each pairing channel and the eigenvectors are 
the pairing form factors. 

The results for 10% electron doping: The five dis- 
joint Fermi surfaces at this doping level are schemati- 
cally illustrated in Fig. (l)(b). For A^ = 32 (the number 
of discretized BZ patches, see Fig. (4)) the few lowest 
eigenvalues of Vf^f as a function of the RG running cut- 
off ln(Ao/A) is shown in Fig. (2)(a). Here Aq, A are the 
initial and the running energy cutoffs. As ln(Ao/A) in- 
creases the most attractive pairing channels has extended 
s-wave symmetry si. In particular, the sign of super- 
conducting order parameter is opposite on the electron 
and hole pockets. The associated gap function, /a(k) 
(here a labels the Fermi surfaces), is shown in Fig. (2)(b), 
and Fig. (l)(a) is a schematic representation of it. Inter- 
estingly, |/a(k)| shows substantial variation on the elec- 
tron pockets (2 and 3). This variation can be under- 
stood qualitatively by considering the following form fac- 
tor Ai + A2{coskx + cos ky). For example, if we expand 
{kx^ ky) around (tt, 0), via kx = Qx^^^^y = Qyi the above 
form factor becomes Ai — 742(cos qx — cos qy). As a result, 
it contains a d-wave like component as q goes around the 
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FIG. 2: (a) The A^ = 32 RG evolution of the scattering am- 
pUtude in the top four attractive Cooper channels for 0.1 elec- 
tron doped system(6.1 electrons per site). Here Ao = 2.2eV 
is the initial energy cutoff, and A is the running energy cutoff. 
The constant a = — l/ln(0.97). (b) The gap function asso- 
ciated with the most attractive pairing channel determined 
from the final (aln(Ao/A) = 101) renormalized interaction 
vertex function. It has a full gap and exhibits significant am- 
plitude modulation on the electron pockets (Fermi surface 2 
and 3). More importantly, the sign of the gap function is 
different on the electron and hole packets. 



electron pockets. If A2 is sufficiently big such form factor 
can exhibit nodes on the electron pocket despite the fact 
it has s-wave symmetry. Although for the parameters we 
have studied the gap function has no nodes on any Fermi 
surfaces, we can not rule out the possibility that for other 
parameter choice the s-wave gap function can have nodes 
on the electron pockets. In Fig. (2) (a) we also show the 
RG flow of the three leading sub-dominant pairing chan- 
nels dx2_y2^ extended-5 with nodes 52, and dxy. For the 
parameter and doping we have studied, triplet pairing is 
not favored. 

Most significantly, from our calculation the pairing 
mechanism can be determined. By monitoring the RG 
evolution of the vertex function, the following sequence 
of events are observed. First of all, the bare Hamiltonian 
contains interaction that can drive SDW. Upon RG these 
scattering vertices grow stronger. Among these grow- 
ing scattering amplitudes there are ones that also cause 
inter-pocket pair scattering as shown in Fig. (l)(b). They 
serve as the "seed" for the growth of other inter-pocket 
pair scattering processes in subsequent RG steps. It is 
important to note that the inter-pocket pair scattering 
drives pairing even when it is positive [21]. In that case 
the superconducting order parameter will have opposite 
sign on the two Fermi pockets, as in our extended s-wave 
pairing case discussed above. From these observations we 
conclude that pairing is driven by the antiferromagnetic 
correlation. 

The results for 10% hole doping: The Fermi sur- 
face topology for this doping value is the same as that 
of 10% electron doping (Fig. (l)(a)). The few lowest 
eigenvalues of the Vff matrices as a function of the RG 
running cutoff is shown in Fig. (3) (a). Again the most 
attractive pairing channel is extended 5-wave like. The 
gap function associated with it is shown in Fig. (3)(b), 
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FIG. 3: (a) The N — 32 RG evolution of the scattering am- 
plitude in the top four attractive Cooper channels for 0.1 hole 
doped system(5.9 electrons per site), where Ao = 2.2ey and 
a = — l/ln(0.97). (b) The gap function associated with the 
most attractive (s-wave) pairing channel determined from the 
final (aln(Ao/A) = 101) renormalized interaction vertex func- 
tion. It has a full gap, and exhibits significant amplitude 
modulation on Fermi surfaces 0, 1, 2 and 3. 



and similarly represented by Fig. (l)(a). Compared to 
Fig. (2)(b) the gap functions on the (tt, tt) hole packets 
is considerably smaller. The mechanism for the electron 
doping applies to this case as well. 

Thus for both electron and hole doped cases discussed 
above extended 5-wave is the pairing symmetry. In con- 
trast, non- s-wave symmetry such as d-wave or triplet 
pairing are not favored. Based on the above results, we 
predict that the extended s-wave is the pairing symmetry 
of the FeAs based superconductors. 

We have also tried several other sets of interaction 
parameters, for instance, varying Jh over the range 
< Jh < ^1,2, varying the ratio and the overall mag- 
nitudes of /7i^2- In addition we have varied the doping 
while maintaining the topology of hole Fermi pockets at 
r and the electron pockets at M. The result that the 
gap function assumes s-wave symmetry and takes oppo- 
site phase on the electron and hole pockets remains un- 
changed. This is even so for significant electron doping 
so that the hole pocket at (tt, tt) disappears. The quali- 
tative nature of our results is not affected by the choice 
of the initial cutoff, as long as it is comparable to the 
band width. Interestingly, the degree of gap variation on 
the electron and hole pockets does depend on interaction 
parameters and doping. For relatively weak interaction 
parameters, we have tried Ui = 2.4eV, U2 = 1.2eV and 
Jh = 0.42eV, the results remains qualitatively similar. 
However, for Ui = 1.2ey, U2 = 0.9eV and Jh = O.lSeV, 
we have not observed any divergence of any scattering 
vertices to the lowest temperature (O.lmeV) we have 
studied. 

FRG details: The essential complication of our FRG 
in comparison with that in Ref. [17] is the fact that the 
Fermi surface(s) (if any) associated with each band is dif- 
ferent. This requires us to discretize the BZ differently for 
different bands (see Fig. (4)). For the band with electron 
pockets, the BZ is chosen to be two 45°-rotated squares. 
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FIG. 4: The Brillouin zone patching scheme. Here the five 
disjoint Fermi surfaces are labeled from to 4. The left panel 
is for the band with smaller hole pocket at (0, 0) (FS 0). The 
right top panel is for the band with larger hole pocket at (0, 0) 
(FS 1) and hole pocket at (tt^tt) (FS 4). The right bottom 
panel is for the electron pockets (FS 2 and FS 3). 



centered at (0, — tt) and (tt, 0) (see the bottom panels of 
Fig. (4)). Each square is divided into N/2 = 16 patches 
by radial lines from its center. The BZ discretization for 
the band containing the large hole pocket at (0, 0) and 
(tt, tt) is similar, but with the two 45°-rotated squares 
centered at (0,0) and (7r,7r). The BZ discretization of 
the band containing the small hole packet at (0, 0) is the 
same as what is used in Ref. [17]. 

Like in Ref. [17] the value of the renormalized ver- 
tex function at a given (ki, k2, ka) in the remaining (un- 
integrated) BZ region is made equal to the value of such 
function with ki,k2,k3 projected to the closest Fermi 
surface segments. As a result, for the three bands with 
five disjoint Fermi surfaces each discretized into M seg- 
ments (in our case M = 16), we need to recompute 
1.536 X 10^ interaction vertices at each step of the RG. 
The top and bottom bands that have no Fermi surface are 
ignored. Due to the existence of many degenerate points 
in the bandstructure, the gauge choice of the numerical 
Bloch wavefunctions needs to be fixed with caution. 

Conclusion We have performed a numerical func- 
tional renormalization group calculation to determine the 
pairing symmetry of the FeAs superconductors. We find 
5-wave with the opposite sign on the electron and hole 
pockets is the pairing symmetry for both electron and 
hole dopings. By monitoring the growth of the renor- 
malized vertex function, we find the pairing mechanism 
is driven by the antiferromagnetic correlations. We hope 
this prediction will be scrutinized by future experiments. 
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